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Abstract 

We study the dynamics of flat Friedmann-Robertson- Walker (FRW) cosmologies in 
the presence of a triplet of massive vector fields with SO (3) global symmetry. We find an 
i£ 3 -symmetric ansatz for the vector fields that is compatible with the i? 3 -invariant FRW 
metric and propose a method to make invariant ansatze for more general cosmological 
models. We use techniques of dynamical systems to study qualitatively the behaviour of 
the model and find, in particular, that the effective equation of state of the system changes 
gradually from a radiation-dominated to a matter-dominated form and that the scale of 
the transition depends on the mass of the gauge fields. 
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I. INTRODUCTION 



Cosmological solutions to the Einstein- Yang-Mills (EYM) equations have been known 
for some time. In particular, solutions corresponding to open, closed and flat Friedmann- 
Robertson- Walker (FRW) universes, in the case where the gauge group is SU(2), were 
studied in refs. [1,2]. The generalization of these results to arbitrary gauge groups, for 
closed FRW models, was done in refs. [3,4] (the euclidian case), where an 5'0(4)-symmetric 
ansatz is derived from the theory of symmetric fields on homogeneous spaces [5]. This 
S'0(4)-symmetric ansatz has been also used for constructing the wave- function of a Uni- 
verse dominated by radiation [6] and to study cosmological implications of generalized 
Kaluza-Klein theories [7]. The case of a flat FRW inflationary model in the presence of 
i±? 3 -symmetric gauge and scalar fields was studied in ref. [8]. 

In treating the classical dynamics of the very early universe with the help of a vari- 
ational principle applied to actions derived from particle physics theories, one hopes to 
obtain a better understanding of the effective equations of state valid in the classical 
regime. On the other hand, this approach may be useful in order to bridge the gap be- 
tween the rather distinct descriptions of the classical and quantum regimes that have been 
used sofar. In this context, it is quite natural from the particle physics viewpoint to study 
non-abelian gauge theories in a cosmological setting. In realistic particle-physics-motivated 
cosmological models, however, the vector bosons normally acquire a mass thereby breaking 
the gauge symmetry. As this process takes place in the very early Universe, it is of interest 
to investigate cosmological solutions of the resulting theory. 

The purpose of this paper is to investigate homogeneous and isotropic solutions to 
the EYM system after symmetry breaking. We shall assume that the gauge symmetry 
one starts with is SOi (3), which subsequently breaks down to a global SOi (3) (through 
a mechanism which we choose not to specify but can be thought of as the usual Higgs 
mechanism) and that, in this process, the gauge fields acquire a mass. We derive an E 3 - 
symmetric ansatz for the triplet of massive vector fields with SOi (3) global symmetry 
that is compatible with the .E 3 -invariant (flat) FRW ansatz for the metric. We compare 
the dynamical behaviour of the system before and after symmetry breaking and find, in 
particular, that, in the latter phase, the equation of state changes with time, gradually 
from a radiation-dominated to a matter-dominated form and that the scale of the transition 
depends on the mass of the vector field. 

Homogeneous and isotropic solutions exist for other gauge groups and symmetry 
breaking patterns (see Appendix). For the case of electroweak symmetry breaking, SU (2)® 
U(l) down to a U(l) however, at least for the ansatz proposed below, homogeneity and 
isotropy cannot be achieved, as the vector fields acquire different masses after symmetry 
breaking thus preventing the possibility for the energy-momentum tensor (EMT) to have 
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the perfect fluid form. 

Recently, there has been great interest in models where inflation is driven by massive 
vector fields [9]. In these works, the vector field is a U(l) gauge field which, upon symmetry 
breaking, acquires a potential leading to inflation. Such a vector field is not allowed 
by an isotropic spacetime and consequently the authors have to choose an anisotropic 
cosmological model (Bianchi I), which, at least in some cases, leads to excessive anisotropy 
at late times. The present work shows that if one starts from a non-abelian gauge theory 
this problem may not arise since, in this case, it is in general possible to built models which 
are compatible with an FRW cosmology from the start. However, the presence of a mass 
term does not by itself lead to sufficient inflation and more involved potentials have to be 
used in order to achieve it [9,10]. In this paper, we shall be interested only in analysing the 
consequences of a mass term for the vector fields and therefore we do not find inflationary 
solutions. 

The paper is organized as follows. In section 2, we derive the £' 3 -invariant ansatz for 
a triplet of massive vector fields with SOi(3) global symmetry. In section 3, we briefly 
discuss non-trivial solutions to the coupled EYM equations with gauge group SOi(3) which 
correspond to the spatially flat FRW universe. We then examine the case where the vector 
fields become massive and derive the corresponding equations of motion from an effective 
action built upon substitution of the i? 3 -symmetric ansatz in the initial action. In section 
4, we analyse the solutions of the equations of motion using methods of the theory of 
dynamical systems and in section 5 we present our conclusions. Finally, in the Appendix, 
we present a method to find G-symmetric ansatze (where G is the isometry group) for 
massive vector fields with an arbitrary global internal symmetry. 

II. ANSATZ FOR VECTOR FIELDS WITH SO(3) GLOBAL INTERNAL 
SYMMETRY 

The spacetime in flat FRW universes has the form 



where the six-dimensional euclidean group E 3 is the isometry group of the spatial hyper- 
surfaces 



The isotropy group SO(3) leaves the origin o E R 3 = E 3 /SO(3) invariant but rotates the 
curves which pass through o and therefore defines a representation of SO (3) in the tangent 
space at the origin T R S . This is the so-called isotropy representation which, in this case, 
is just the vector representation 3 of SO(3). 



M 4 = R 4 = Rx E 3 /SO{3), 



(1) 



M 3 = R 3 = E 3 /SO(3). 



(2) 
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The fact that there are no singlets in the isotropy representation has far reaching 
implications. Indeed, if a vector field X in R 3 is i? 3 -invariant then, in particular, it must 
be invariant at o, i.e., a singlet under the action of SO (3) on T Q R 3 . Since this is not 
possible, we conclude that there are no i? 3 -invariant vector fields in R 3 . 

Being physically obvious for R 3 (as it just means that there are no homogeneous and 
isotropic vector fields) this result is important as it applies for any homogeneous space G/H 
for which the isotropy representation has no trivial singlets. This is also the case for the 
spatial hypersurfaces S 3 and H 3 of closed and open FRW models which, therefore, do not 
have vector fields invariant under the action of the corresponding isometry group either. 
This explains why a single vector field is incompatible with FRW geometry; indeed, being 
non-spatially homogeneous and isotropic, the vector fields generate anisotropic energy- 
momentum tensors which therefore cannot be proportional to the Einstein tensor of the 
FRW metric. 

In this section, we shall study the case of flat FRW models, leaving the analysis 
of the more general case of a (not necessarily four-dimensional) spacetime of the form 
M d = R x G/H to the Appendix. Models with d > 4 are relevant to the study of the 
cosmological implications of Kaluza-Klein theories [7] . 

If, instead of a single vector field, we consider a theory with a multiplet of (covariant) 
vector fields A a ^a = 1, ...,N, where a is an internal index, we may hope that, although 
the EMT corresponding to each field A a ^ is not £' 3 -invariant, it would be possible to relate 
the asymmetries of the different vector fields in such a way that the total EMT is E 3 - 
invariant. Let us then consider the case of a triplet of massive vector fields with global 
SO i (3) symmetry. The generators of the isometry group G = E 3 satisfy the following 
commutation relations 



[Ti,T j }=0, 

\Qii Qj] = e ijkQki (3) 
[Qii Tj] = CijkTk', 2, J, k = 1, 2, 3 

so that the Ti span an (abelian) invariant subalgebra of Lie(i? 3 ) — the algebra of transla- 
tions. The corresponding Killing vector fields are 



where x % are the cartesian coordinates in R 3 and -J^ together with form a (global) 
moving frame in R 4 = Rx E 3 / SO(3), {X^} — {J^-, -£~r}, with the dual coframe being given 
by {uj^} = {dt,dx 1 } {p, = 0, 1, 2,3). The remaining Killing vector fields, corresponding to 
the rotations Qi, read 
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d 

dx k ' 



(5) 



so that 



C x .Xj = [Xi,Xj] = 0, (6 -a) 

A^- = -C Xj Yi = \Yi, X 3 ] = e ijk X k , (6 - 6) 

LyYi = [Y h Yj\ = e ijk Y k . (6-c) 

The conditions of spatial homogeneity and isotropy read, infmitesimaly 



C Xi A = 0, (7 -a) 

C Yi A = 0, (7-6) 

where A = A®(t, x)uj^L a , the L a being the generators of the internal group 50/(3) 

[L a , Lb] = e abc L c . (8) 

In agreement with our previous considerations, we conclude that the only solution of (7) is 
the trivial one, i.e. such that only the zeroth component of the covector fields is different 
from zero: A = A®(t)dtL a . Condition (7-b) is, however, too restrictive and we now consider 
a weaker version of it by allowing the field A to be such that infinitesimal spatial rotations 
generate internal SOi(S) rotations 



C x ,A = 0, (9 -a) 

C Yi A = -[L i ,A}. (9-6) 

It is clear that, if (9) takes place, then, although A is not i? 3 -invariant, it generates an 
EMT that is invariant under the action of E 3 . Indeed, the action of E 3 in the EMT is 
equivalent to an internal rotation with respect to which the EMT is invariant 

C Yz EMT(A) = -5 Lz EMT(A) = 0. (10) 

The conditions for spatial homogeneity, eq. (9-a), imply that the components A a ^ depend 
only on the time variable t. In order to solve (9-b), let us consider the linear space of vector 
fields with constant coefficients: C = {X = a M X M ,a^ G R}. Eq. (6-b) implies that X e C 
transforms as 1 ©3 under SO (3) spatial rotations. On the other hand, the one- form A can 
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be considered, for any t, as a linear mapping from C to the internal space (i.e. to the Lie 
algebra of 50/(3)) 



A(X) = A(a»Xj = A«(t)a»L a , (11) 
in which the triplet representation 3 of 50/(3) is realized. Rewritting (9-b) in the form 

A(C Yi X) = [L i ,A(X)}, (12) 

we see that A intertwines the representations 1 ©3 and 3 of 50(3) and 50/(3) ls ^ m 50(3), 
respectively. Schur's lemma then implies that A vanishes between subspaces of non- 
equivalent representations 

A = <=► A Q = 0, (13) 

and is proportional to the identity operator between subspaces with equivalent represen- 
tations 

A (^~) = Xo(t)Ll ^ A ' {t) = X ° m ' (14) 

where Xo(t) is an arbitrary function. Alternatively, the i? 3 - symmetric ansatz for A satis- 
fying (9) and therefore generating S 3 -invariant EMT, can be written as 

A = Xo (t)dx i L i . (15) 



III. EFFECTIVE ACTION AND FIELD EQUATIONS 

A. Massless case 

We start with the following action for the coupled EYM system with 50/ (3) gauge 
group 

s= L dlxV ^{^ R+ ^ TriF ^ ] )- (i6) 

where k 2 = 8tvG and e is the gauge coupling. Our conventions and notation for the 
gravitational part of the action correspond to those of refs. [3,4]. The gauge field strenght 
is = d^Ajy — d u A^ + [A^A^], where A^ = A a ^L a . We assume that the spacetime 
manifold M 4 is topological^ of the type (1), where {t D } x £ 3 /50(3) = {t } x R 3 , t e R, 
are the spatial orbits of the isometry group G = E 3 = R 3 ® 50(3). 
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The FRW ansatz, for flat models, gives the most general form of an i? 3 -invariant 
metric 



ds 2 = -N 2 (t)dt 2 + a 2 (t) ((dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ) (17) 

where the lapse function N(t) and the scale factor a(t) are arbitrary non-vanishing func- 
tions. 

Cosmological models in the presence of gauge fields have been investigated in some 
detail in the literature [1-3,7,8]; for these models, the local symmetry of the action allows 
the use of the theory of symmetric gauge fields [5] in order to make the necessary ansatze. 
The G-symmetric ansatz corresponds to fields which, under the action of the isometry 
group G, undergo a local internal rotation; however, in theories with global symmetry 
such fields do not, in general, lead to G-invariant EMTs, as can be seen from eq. (10) (for 
more details, the reader is referred to the Appendix). Ansatz (15), for (covariant) vector 
fields with global SOj(3) internal symmetry is, of course, also valid in the case of action 
(16), with local SOi (3) symmetry, and we shall use it in what follows. It is easy to check 
that this ansatz leads to a (traceless) energy-momentum tensor (EMT) with the required 
E 3 -invariance 



T ° = -p(t), TS=p(t)8i j , 



(18) 



i.e., it corresponds to the EMT of a perfect comoving fluid with energy density p and 
pressure p given by 



V = 



Vgf 



Aa 2 e 2 N 
1 

Aa 2 e 2 N 
1 



Y 2 iV ! 



Y 2 iV ! 



T 2 

r 2 

>2~ 



(19) 
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where the dot denotes the derivative with respect to t; moreover, the equation of state for 
the system is clearly that of a radiation fluid: p = 3p. 

There are, for this model, two equivalent methods to obtain the dynamical equations 
for the functions a(t) and Xo(t): one can either substitute the ansatze, eqs. (13), (14) and 
(17), into the field equations derived from action (16) or substitute the ansatze directly 
into the action and then get the dynamical equations from this effective action. The latter 
method has the advantage of simplifying the study of the equations of motion and we shall 
adopt it in what follows. The effective action is given by (disregarding the infinite volume 
of the spatial hypersurface) 
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/•*2 

S ef f = dt 



3 d 2 a 3a f xl N 



k 2 N + 4e 2 N \ 2 a 2 Vgf 



(20) 



The equations of motion are obtained upon variation of this action with respect to N(t), 
a(t), and x (t). We get, respectively, in the "gauge" N = 1 



h 2 
A; 2 

H + H 2 = - — (p + 3p) 



Xo + XoH 



6 

'2a 2Xo ' 



(21 - a) 
(21 - b) 
(21 -c) 



with p and p given by eq. (19) and H = a/a. 

The Yang-Mills equation, (21-c), takes a much simpler form in the conformal time 

[3,4] 



dr] 



dt 



(22) 



since, in this case, this equation decouples from the remaining ones. Indeed, substituting 
(22) in (21-c), we obtain 



II _ _Xo 

Xo - 2 , 



(23) 



where the prime denotes the derivative with respect to the conformal time. The solution 
of eq. (23) is given by 



Xo = f(v,E,r} ), 
where / = f(rj; E, t]q) is a function defined implicitly by 

V -7 ]0 = (2E)- 1/4 F( 1 (f,E),r). 
In the above equation, F(x, k) denotes the elliptic function of the first kind [11], 

7 (/,i?) = cos- 1 [(8£)- 1/4 / , 
V2 

and the constant E is the conserved "mechanical energy" of the solution 



(24) 



(25) 



(26) 
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E=l( X 'o) 2 + V gf . 



(27) 



Using a "mechanical" analogy, we can consider eq. (23) as describing the motion of a 

'2 

particle with kinetic energy 4*-, potential energy V g f and total energy E. The above 
solution describes therefore a particle oscillating in a potential well between the turning 
points: -(8E) 1 / 4 and (8E) 1 / 4 . 

B. Massive case 



Consider now the action which results from adding a mass term for the gauge fields 
[12] to the action of eq. (16) 

S' = d^x^-g (±R + ^Tr [F^F^] + \ m 2 Tr[A,A^ . (28) 

Such a mass term would arise, e.g., via the Higgs mechanism, with a complex doublet 
of scalar particles acquiring a vacuum expectation value. The above action is no longer 
gauge invariant and only a global SOi(3) symmetry remains. Substituting the ansatze, 
eqs. (13), (14) and (17), into the energy-momentum tensor corresponding to the action 
(28), we obtain 

T ° = ~P- \m* 
T t J = (p-^m 2 

with p and p given by eq. (19), which shows that the perfect fluid form of the EMT is 
preserved in spite of the presence of the mass term. Notice that, if instead of a mass term 
we had a quartic potential, this could be absorbed in a redefinition of V g f, which would 
lead us back to the results of the previous subsection. 

We are interested in spatially homogeneous and isotropic cosmologies associated with 
the action (28). Hence, we restrict ourselves to the set of ^-symmetric configurations 
for the fields g^ v and A^. Substituting these ansatze into (28), we obtain for the effective 
action 



Xo 

a 



— =-P. 



(29) 



S, 



eff 



\ dt 



3 a a 



+ 



3a 



k 2 N 4e 2 N \ 2 



N 2 



a" 



^ " -^V gf - -Nm\^a 



(30) 



The equations of motion, derived upon variation with respect to N(t), a(t) and x (t) are, 
respectively, in the "gauge" N = 1: 
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H 2 = y p, (31 - a) 

H + H 2 = -^(p + 3p), (31-6) 
Xo + XoH = -^Xo-2m 2 e 2 Xo , (31 - c) 

with p and p as given by eq. (29). Consider now the following change of variables 

r,Xo 

X — K , 

a 

- fc2 *° 

V ~ 2V2 a ' (32) 
z = kH, 
t = t/k. 

The system of eqs. (31), in the new (dimensionless) variables, can be written as 



x T = 2\p2y — xz, (33 — a) 

Vt = -^-p 2 x - 2yz - -^£ 3 , (33-6) 

z T = -2z 2 + j^V, (33 - c) 

z 2 = y 2 + ^x 4 + ^p 2 x 2 , (33-d) 

where p 2 = m 2 k 2 and we have set e = 1; the index r denotes the derivative with respect 
to r. Eqs. (33-a)-(33-c) define a three-dimensional dynamical system in the variables x, 
y, z and eq. (33-d) is a constraint equation which defines, in i? 3 , the phase space of the 
dynamical system. Substituting the constraint into eqs. (33-a), (33-b), we are left with a 
two-dimensional dynamical system in the variables x, y 

x T = 2V2y -x(y 2 + —x 4 + -p?x 2 J , 

V J 1 (34) 

\/2 9 n ( 5 , , 1 4 , 1 2 2^ 2 v^2 



Vt = — x ~ 2y y + 32 x + 4^ x J ~ ~Y X ' 

where we have considered the positive root of the constraint (33-d) as we are interested in 
expanding models. 
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It is this system of differential equations (34) that we are going to study using the 
methods of the qualitative theory of dynamical systems [13]. 



IV. QUALITATIVE ANALYSIS 

The dynamical system (34) has just one critical point in the finite region of variation 
of x, y, the origin, hereon referred to as F(0, 0) or simply F. The remaining critical points, 
if any, will lie at infinity. In order to study them, we complete the phase space with an 
infinitely distant boundary (for which x 2 + y 2 = +oo). It is also convenient to perform 
a change of variables to polar coordinates, (x, y) — > (r, 0), followed by a compactification 
of the entire phase space into a circle of unit radius. This can be achieved through the 
introduction of yet new radial and time coordinates p and (, such that 

(0<P<1); d i = rA^- (35) 



l-p " " dr (1-p) 

In the variables (p, 0, £), the system (34) becomes 

u 2 

p c = V2(2 - ^-)p(l - p) 3 sin 9 cos 9 - p 2 (1 - p) (1 - sin 2 9) f(p, 9) 

- ^V(l - p) sin0cos 3 9 = U(p, 9), 
$ c = -y/2(l - p) 2 [2 sin 2 9 + ^ cos 2 O^j - ^-p 2 cos 4 9 - psin9 cos f(p, 9) 

where f 2 (p, 9) = (1 - p) 2 sin 2 9 + fj cos 4 9 + ^(1 - p) 2 cos 2 9. 

The critical points in the infinitely distant boundary (p = 1) are the solutions of the 
equations: 11(1,0) = and #(1,0) = 0, namely, iVi(l,§), iV 2 (l,^), Si(l,^), S 2 (l, I f) 
— see Fig.l. 

We now turn to the analysis of the nature of the critical points of the system. First, 
we would like to point out that, since (36) is invariant under — > + tt, we need only 
consider two of the four critical points at infinity; hence, we are left with three points to 
analyse, F(0,0) and, e.g., N^l, f) and Si(l, *f). 

Regarding F(0,0), it is easy to check that this point is degenerate (see Table 1) and 
therefore the Hartman-Grobman theorem [13] does not apply. We resort to the method 
known as "blow-up" [13], which allows us to establish that F is a focus. In order to 
investigate its stability, we use Liapounov's theorem [13]. Consider the Liapounov function 

£(x,y)=2(y 2 + ±x* + ^x 2 ). (37) 
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Since C(x, y) obeys the conditions 



£(0,0) = 0, £((x,y)^(0,0))>0, 
<0 in R 2 /{(0,0)}, 



dC _ . n2ir/n _ (38) 



dC 

we conclude that F is a stable focus and all trajectories in the x,y plane asymptotically 
approach this point. Furthermore, trajectories in the neighbourhood of F are clockwise 
directed spirals. 

Regarding the critical point at infinity Si, the linear approximation of the system (36) 
around this point is sufficient and we conclude that Si is a saddle point (cf. Table 1). As 
for N 1: the situation is similar after we perform a change of time variable (( — > () such 
that ^ = (p — 1) and we conclude that N\ is an improper unstable node (cf. Table 1). 
We shall also use the fact that the trajectories in the vicinity of Ni are tangent to the 
eigenvector which corresponds to the lowest eigenvalue — see Table 1 and Fig. 1. 

In order to get a more detailed picture of phase space, we divide it into four distinct 
regions and find approximate analytic expressions for the trajectories in each of these 
regions, namely 

^ 2 1 
a x = — cos r 

' fir + d 



1 

b) x + y = 2C 2 x 2 I x | ~| y |> 1 

c) y = x + C 3 | x | >| y |> 1 

„2 I 1, I nj. „2 



(39) 



d) y = C A x 2 | V | > x z ; | x |, | y |> l 

where the Ci (i = l, ...,4) are arbitrary constants. In case a), i.e., near the origin, we see 
that the shape of the trajectories depends on the value of \i 2 . In fact, for small values of 
/j, 2 the trajectories wind around F with an oblate shape whereas for large values of \i 2 they 
have a prolate shape. It is indeed in this region that the behaviour of the system is most 
sensitive to the presence of the mass term. 

We now turn our attention to the analysis of the effective equation of state for the sys- 
tem near the critical points. In the vicinity of Si, 6 ~ introducing this approximation 
into the dynamical system, we get the following solutions 

On 

a = a Q t 1 / 2 , Xo = -—^j- 4 (f-0) (40) 

(2k ) 

Eq. (29) then implies p ~ 3p, indicating that the effective equation of state near Si is 
approximately that of a radiation fluid. In fact, this is true everywhere except in the 
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vicinity of F. This is to be expected since, were it not for the mass term, the equation 
of state would be p = 3p throughout all regions of phase space and, as we have already 
mentioned, it is near the origin that the mass term most affects the trajectories. 

The effective equation of state near the focus can be analysed using a procedure similar 
to the one explained above for Si , which leads us to conclude that 

a~t 2/3 , Xo~t" 1/3 sin-^= (*->+oo), (41) 

implying that p ~ 0, corresponding to the equation of state for dust matter. 

In addition, we find that the time of transition between the regimes described by the 
different equations of state depends on the value of m 2 such that the transiton occurs earlier 
for more massive vector fields. Before the transition the energy density and pressure evolve 
as expected in a radiation-dominated Universe and, as the system approaches the focus F, 
the Universe becomes matter-dominated and the pressure starts to oscillate around zero 
(see Fig. 2) with ever decreasing amplitude and quickly vanishes. 



5. CONCLUSIONS 

We have shown that the Einstein- Yang-Mills theory, with gauge group SOi(S), broken 
by the presence of a mass term for the gauge fields in such a way that only a global £0/(3) 
symmetry remains, admits homogeneous and isotropic flat FRW solutions. To prove this, 
we have constructed an £' 3 -symmetric ansatz for the (massive) vector fields, parametrized 
in terms of an arbitrary function of time Xo(t)- 

We have studied the equations of motion of this theory using the methods of the 
theory of dynamical systems and found that the system has just one stable critical point, 
the origin, which is a focus. We have checked that the behaviour of the system is basically 
the same as for the unbroken phase except near the origin; in this region, the change is 
such that the effective equation of state of the system becomes that of a matter-dominated 
universe 1 . 

APPENDIX 

Let us consider a spacetime of the form 

1 When the present paper was already in press we became aware of ref. [14] where an 
exact cosmological solution for an isotropic gas of collisionless massive particles has been 
found. This solution also describes a gradual change of the equation of state from p = l/3p 
to p = 0. 
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M d = Rx G/H, (A.l) 

where {t Q } x G/H are spatial hypersurfaces, G is the isometry group and t G R is a time- 
like coordinate. There are two classes of models for which vector fields can be compatible 
with the above geometry, in the sense that solutions of the coupled Einstein-vector field 
equations can be found in a rather simple way: non-abelian gauge theories and theories 
with a global internal symmetry. In the former case, since the theory is invariant under 
an internal local symmetry, the total EMT is also gauge invariant, which implies that, for 
the total EMT to be G-invariant, the fields A® need only to be G-invariant up to a gauge 
transformation. For a gauge group K containing an H subgroup such fields often exist, 
which has allowed the study of FRW cosmologies in the presence of non-abelian gauge 
sectors [3,7,8]. The latter case has, however, to our knowledge, not yet been studied in the 
literature. In section 2, we analysed a particular example, with G/H = E s /SO(3) and 
K = SO i (3). We now study the general case. 

Consider configurations A^(x) for which the action of g G G can be compensated by 
an internal global K-rotation 



g* {A«(x)dx») = ^(Xig-^tA^dx^ (A2)) 

where A : G — > K is an homomorphism and W is the representation of K in the internal 
space. Clearly, configurations satisfying (A. 2) will generate G-invariant (spatially homoge- 
neous and isotropic in the FRW case) EMTs. Notice that condition (A. 2) is less restrictive 
than the one of G-invariance (corresponding in (A. 2) to X(g) = e, 6 G) but still much 
more restrictive than the one for gauge theories since the compensating internal rotation, 
ty(\(g)), cannot depend on the points x G R x G/H. Nevertheless, we will show that, for 
a class of homogeneous spaces G/H which includes the spatial hypersurfaces of closed and 
flat FRW models, field configurations satisfying (A. 2) can be found. 

Let {T~}^"^ G be a basis in Lie(G), the Lie algebra of the isometry group G, satisfying: 



i' J 



= C~T~ and \X~\ be the associated Killing vector fields 
ij k 1 lJ 6 



d 

X^(x) = — (exp(-e T-)x) | e=0 . (A. 3) 

We shall assume that the homogeneous space G/H is such that there exists in Lie(G) an 
invariant subalgebra (ideal) with basis {Ti}f^ G ^ H 

[T i ,T j ]=C^ k , (AA) 

and such that {Xi}f^ G ^ H is a (local) moving frame in G/H. This is obviously the 
case for flat (see sect. 2) and closed FRW models. Consider now the moving frame in R x 
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G/H, {X^} d ^™ Q G/H = {^|, X i }'f™ G ' H , and the moving coframe {uj^} d ^™ Q G/H = {dt, uj 1 ; i = 

1, dimG/H} dual to {X^} d ^™ G ^ H . The commutation relations define a representation 
of Lie(G) acting in the subalgebra of vector fields 

C = {X = c« M , e R}, (A5) 

according to 



where 



Cx.X = a'^X^ {A.6) 



a'° = 0, a' 1 = Cjja*. (A7) 



In the coframe {w' 1 }, we choose the fields A a ^ to depend only on the coordinate t 

A a = A£(t)uA (A.8) 

As the w M are not G-invariant, the fields A a are not G-invariant either. Infinitesimally, 
condition (A. 2) applied to (A.8) gives 

Cx i A a = -^(X(T i )) a b A b , (A.9) 

or equivalent ly 

A a (£ Xl x)=*(\(T l )) a b A b (X). (A.10) 

These are linear algebraic constraints on A®(t), which imply that, as a mapping from C 
to the internal space (for fixed t) 

A^ : a»^B a = A a ^, (All) 

A a ^ intertwines the representation of Lie(G) in C with the representation of 
A(Lie(G)) in the internal space. From Schur's lemma, we then conclude that there exist 
nontrivial solutions of (A. 10) if and only if there are equivalent irreducible representations 
in C and in the internal space. Therefore, in order to make a G-symmetric ansatz for the 
fields A a , one must: 

(i) Check whether (A. 4) takes place for the given G and G/H and choose the associated 
coframe {cu M }. 

(ii) Choose an homomorphism A of the isometry group G to the internal symmetry group 
K such that there are equivalent irreducible representations in the representations of 
G in C (see (A.6), (A. 7)) and of A(G) in the internal space. 
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(iii) Find the intertwining operators A a ^ between the above two representations and sub- 
stitute them in (A. 8). 
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Point 


Eigenvalues 


Eigenvectors 


Classification 


F(0,0) 






Focus 




V2 
16 

V2 
16 


d 
dp 

a 

89 


Saddle 


iVi(l,f) 


2 

1 


— - 2\/2 — 
dp zvz d9 

a 

89 


Unstable node 



Table 1. Classification of critical points from the analysis of the dynamical system of eqs. (34) 
and (36). 
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FIGURE CAPTIONS: 



Fig. 1 All trajectories of the phase diagram of the model asymptotically approach the focus 
F; Ni, A^2, Si, S2 are critical points at infinity. Near the focus, trajectories change 
from oblate to prolate shape as the mass of the vector bosons increases. 

Fig. 2 3p/p as a function of time showing the change of the equation of state of the Universe 
from radiation-dominated to matter-dominated. 
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